ALGEBRAIC-GEOMETRIC CHARACTERIZATION OF AUTONOMY 
AND CONTROLLABILITY OF BEHAVIOURS OF SPATIALLY 

INVARIANT SYSTEMS 



AMOL SASANE 

Abstract. We give algebraic-geometric characterizations of the properties of autonomy 
and of controllability of behaviours of spatially invariant dynamical systems, consisting of 
distributional solutions w, that are periodic in the spatial variables, to a partial differential 
equation 

f_d_ ... J_ d\ w = Q 

\dxi ' ' dx d ' dt J 
corresponding to a polynomial p £ C[£i, . . . , £<j, r]. 



1. Introduction 

Consider a homogeneous, linear, constant coefficient partial differential equation, in 
described by a polynomial p £ C[£i, . . . , r]: 

/ a d d\ 



9 d d 

P ' 



(That is, the operator 

9xi ' ' dxd dt y 

is obtained from the polynomial p £ C[£i, . . . , £d, r] by making the replacements 

6 7T- for k = 1, . . . , d, and r ~> ^. 

In the behavioural approach to control theory, the "behaviour" 53yy;(j>) associated with p 
in W (where W is an appropriate solution space, for example smooth functions C°°(lR rf+1 ) 
or distribution spaces like V'(K. d+1 ) or 5'(lR d+1 ) and so on), is defined to be the set of all 
solutions w G W that satisfy the above PDE (II. lj) . and one characterizes algebraically (in 
terms of algebraic properties of the polynomial p) certain analytical properties of 25yy(p) (f° r 
example, the control theoretic properties of autonomy, controllability, stability, and so on). 
See for example [8], [9], pQ for distinct takes on this in the context of systems described by 
partial differential equations. 

For example let us consider the property of "autonomy", which means the following. 
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Definition 1.1. Let W be a subspace of T>'(M. d+1 ) which is invariant under differentiation, 
that is, for all w G W, 

d d 

w G W for k = 1, . . . , d, and — w G W. 



If p G C[£i, . . . , (la, r], then the behaviour (of p in W) is 

We call the behaviour s 3y\>(p) autonomous (with respect to W) if the only w G Q3yy;(p) satis- 
fying w\t<:o = is io = 0. 

The following Theorem 3.4] is an easy consequence of p. 310, Theorem 8.6.7]. 

Proposition 1.2. Le£ p G C[£i, • •• ,£<i, t] 6e a nonzero polynomial. Then the behaviour 

w : = {» 6 : » (sr - • • ■ • £~A) » - °} 

corresponding to p is autonomous if and only if degp = degp(0, r). 

Here by deg(-), we mean the total degree, which is the maximum (over the monomials 
occuring in the polynomial) of the sum of the degrees of the exponents of indeterminates in 
the monomial. Also, p(0, r) denotes the polynomial in C[r] obtained from p G C[£i, • • • , £d, t] 
by making the substitutions i— >■ for k = l,...,d. 

There has been recent interest in "spatially invariant systems", see for example [3], [3], 
where one considers solutions to PDEs that are periodic along the spatial direction. So it is a 
natural question to ask what the analogue of Proposition II. 21 is. when we replace the solution 
space V(M. d+1 ) with one that consists only of those solutions that are periodic in the spatial 
directions. In this article, our first main result is the following one, characterizing autonomy 
of spatially invariant systems. 

Theorem 1.3. Let p G C[£i, . . . ,(,d, T ] be a nonzero polynomial and A := {ai, . . . ,ad} be a 

linearly independent set vectors in M. d . Then the behaviour 

s^+,(ri == {» e :r(ir l .--.£- d 4) w = 

is autonomous if and only if (V(C^(p))) f](2iriA~ l Z d ) = 0, where A is the matrix with its 
rows equal to the transposes of the column vectors ai, . . . , a^: 

T 



.4 



a d 

Here T>{(R d+1 ) is, roughly speaking, the set of all distributions on R d+1 that are periodic 
in the spatial direction with a discrete set A of periods. The precise definition of T>' A (R d+1 ) 
is given in Subsection ll.il Also, in the above result, the condition 

(V(Ct(p)))(~)(27riA- 1 Z d ) = ® 

is an algebraic-geometric condition, saying that the variety V(C^(p)) of the "^-content" of 
p does not meet the discrete set of points in 2iriA~ lr L d . The notion of the ^-content of a 
polynomial is defined in Subsection 11.21 
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1.1. The space X>^(R d+1 ). For a G R d , the translation operation S a on distributions in 
V'(R d ) is denned by 

(S a (T), = (T, y>(- + a)} for all 92 G P(R d ). 

A distribution T G Z>'(R d ) is said to be periodic with a period a € R rf if T = S a (T). 

Let A := {ai, . . . , a^} be a linearly independent set vectors in R d . We define T>' A (R d ) to be 
the set of all distributions T that satisfy 

S ak (T)=T, fc = l,...,d. 

From [5j §34], T is a tempered distribution, and from the above it follows by taking Fourier 
transforms that (1 — e 2?nak ' y )T = for k = 1, . . . , d. It can be seen that 

f = Yl «v(T)5 v , 

for some scalars a v G C, and where A is the matrix with its rows equal to the transposes of 
the column vectors ai , . . . , ad : 



A :-- 



a d T 



Also, in the above, 5 V denotes the usual Dirac measure with support in v: 

(5 V ,V) = V(v) for V G P'(R d ). 

By the Schwartz Kernel Theorem (see for instance [3 p. 128, Theorem 5.2.1]), T>'(R d+1 ) is 
isomorphic as a topological space to £(V{R),V'{R d )), the space of all continuous linear maps 
from T>(R) to T>'(R d ), thought of as vector- valued distributions. For preliminaries on vector- 
valued distributions, we refer the reader to [2]. We indicate this isomorphism by putting an 
arrow on top of elements of V'(R d+1 ). Thus for w G V'(R d+1 ), we set w G C(V(R),V'(R d )) 
to be the vector valued distribution defined by 

(w(ip),tp) = (w,i/)® if) 

for (p G V(R) and if) G V(R d ). We define 

V' A (R d+1 ) = {we V'(R d+1 ) : for all ip G £>(R), tZf(^) G ^(R*)}. 

Then for w G V A (R d+1 ), 

-to G V A (R d+1 ) for jfe = 1, . . . , d, and — w G V A (R d+1 ). 



dx k AV ' — df 

Also, for w G X>^(]R d+1 ), we define u; G £>'(M d+1 ) by 

for 99 G 2?(R) and ^ £ X , (R ci ). That this specifies a well-defined distribution in V(R d+1 ), can 
be seen using the fact that for every <3> G D(R d+1 ), there exists a sequence of functions (^ n ) n 
that are finite sums of direct products of test functions, that is, = XlfcV'fc ® ^Pki where 
ip k G V(R d ) and 99*. G T>(R), such that * n converges to $ in £>(R d+1 ). We also have 

w = 2myj c w for k = 1, . . . , d, and ^-w = ■J-'S. 



<9t dt 

Here y = (yi, . . . , yd) is the Fourier transform variable. 
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1.2. The ^-content of a polynomial and its variety. If p is a polynomial in C[£i, . . . , r], 
then we can write 

p = c + cir + c 2 r 2 + • • • + c n r n G C[fi, . . . ,&][r], 

where Co, ci, C2, . . . , c n € C[£i,...,£j. The ^-content C^(p) is defined as the ideal in the 
polynomial ring C[£i, . . . , £d] generated by c , ci, c 2 , . . . , c n . 

Given a set L of polynomials from C[£i, • ■ • , £<j], the notation V(I) means the (affine) variety 
(in C d ) of /, that is, the set of all common zeros of each polynomial in /: 

V(J) = {v G C d : p(v) = for all p € /}. 

Thus given a p = cq + C]T + c 2 t 2 + • • • + c n r n G C[£i, . . . , £d][r], we have 

V(Cf(p)) = {v G C d : c (v) = • • • = c n (v) = 0}. 

1.3. Controllability. 

Definition 1.4. Let W be a subspace of D'(M. d+1 ) which is invariant under differentiation 
and p G C[£i, . . . , £d, t]. Then we call *Byv(p) controllable if for every w;i,w; 2 G <Byv(p), there 
isaw£ 23>v(p) and a r > such that to|/ = ^i|(-oo,o) an d H(t,+oo) = u, 2|( r ,+oo)- 

The second main result in this article is the following one, characterizing controllability of 
spatially invariant systems: 

Theorem 1.5. Let p G C[£i, . . . , £d,r] 6e a nonzero polynomial and A := {ai, . . . , a,}} 6e a 

linearly independent set vectors in M. d . Let 

= {■» € ^ (Ri+1) : p (sr - • • • ■ h It) w = °} * 101 

Then the following are equivalent: 

(1) 2$x>/ (R<j+i) (p) ^ s controllable. 

(2) co(27riv) = • • • = c n (27uv) = /or a// v G jd -1 ^, where for k = 0,1, ... ,n, the 
polynomial cu G C[£i, . . . , £d] is the coefficient of r k in the decomposition of p as 
p = Co + c\T + • • • + c n T n , and A is the matrix with its rows equal to the transposes of 
the column vectors ai , . . . , ad . 

2. Proof of Theorem 11.31 

Before we prove our main result, we illustrate the basic idea behind 'If part: by taking 
Fourier transform, the partial derivatives with respect to the spatial variables are converted 
into the polynomial coefficients afc(27riy), where y is the vector of Fourier transform variables 
y . . . , But the support of w is carried on a family of lines, indexed by n G Z d , in 
parallel to the time axis. So we obtain a family of ODEs, parameterized by n G Z d , and by 
"freezing" an n G Z d , we get an ODE, where for a solution we can indeed say that zero past 
implies a zero future, and so the proof can be completed easily. 

Proof of Theorem ] 1.31 'Only if part: Let v G C d be such that 

v G V(C f (p))P| (^mA^lA . 

Consider w := e v X <S>0, where O is any nonzero function in C°°(K) which has a zero past. Here 
the notation • is used for the usual Euclidean inner product in the complex vector space C^: 
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thus we have vx := v\X\ + - ■ -A-v^Xd, where v = (v\, . . . ,Vd) G C d and x = (x±, . . . ,Xd) G 
Take for example given by 



e(t) = 

Then we have that w G V' A (R d+1 ), since 



if t < 0, 
e" 1 ^ if t > 0. 



S ak «; = e v -( x + a ><) (gi 6 = e vak e vx ® 6 = 1 • e v ' x <g> 9 = tw. 

Moreover, w has zero past, that is, w\ t< o = because Q\t<o = 0- If P = c o + c i r + ' ' ' + c n T n , 
where c ,ci,...,c„ G C[£i, . . . , then c ,...,c n G and so c (v) = ••• = c„(v) = 0. 

Using this, we can see that w G s BD' A (m. d + 1 )(p) because 

* (£• ' ' ' • h I) w = ( Q)(v)e + Cl(v)e ' + ' ' ' + c " (v)e '"') eVX = °' 

Consequently, ©x>' (R d + 1 ) (p) i s n °t autonomous. This completes the proof of the 'only if part. 



Tf part: Suppose that V(Q(p)) f| {2ixiA~ Xr L d ) = 0. Let w G £>^(R d+1 ) satisfy w|t <0 = and 

Let p = c + cit + • • • + c n T n G C[£i, . . . , £d][r], where c , c 1 ,c 2 ,...,c n G C[£i, . . . and 
c n ^ as a polynomial in <C[£i, . . . Upon taking Fourier transformation on both sides of 
the equation (|2.ip with respect to the spatial variables, we obtain 

d / d \ n 

co(27riy)w + ci(2vriy) — w + ■ ■ ■ + c n {2my)^ — J to = 0. (2.2) 

For each fixed (p G £>(M), w G T>' A (R d ), and so it follows that 

w = ^2 oi v (w, tp)5 v , (2.3) 

for appropriate scalars a v (T, 99) G C. In particular, we see that the support of w is contained 
in A~ lr L d x [0, +00). Thus each of the half lines in A~ lr L d x [0, +00) carries a solution of the 
differential equation (|2.2p . and w is a sum of these. We will show that each of these summands 
is zero. It follows from (|2.3p that the map ip 1— > a v (w,(p) : £>(R) — > C defines a distribution 
T v in D'(R). Moreover, the support of T v is contained in [0, +00). From (|2.2p . we see that in 
a small enough neighbourhood of v G A~ lr L d , we have 

co(27rzv)T v + ci(27riv)— T v + • • • + c„(2vriv) (-) T v = 0. 

Owing to our hypothesis that V(C^(p)) (2i\:iA'~ 1 'L d } = 0, we know that at least one of the 
coefficients co(27riv), . . . , c n (2iriv) is nonzero. Together with the fact that ^V^^o) = 0, this 
yields that T v = on R. Hence w = in a nighbourhood of v G A~^7L d . As the choice of 
v G A~ l Z d was arbitrary, (|2.3p implies that w = and so u> = 0. Consequently, the behaviour 
*B-p/^ K d+i)(p) is autonomous. □ 
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3. Proof of Theorem 11.51 

Proof of Theorem HOI (1) (2): Suppose that <Bw(p) is controllable. Let v G A" 1 !/ 1 be 
such that Cfc(27riv) ^ for some k G {0, 1, . . . , n}, and q(27tzv) = for all £ > k. Let r > 0, 
and let 6 G C°°(R) be any solution to the ODE 

co(2vriv)e + ci(27rrv)e' + • • • + c k (2iriv)Q {k) = 0. 

and e| (T)+oo) ^ 0. Define w 2 G £>'(R d+1 ) by w 2 := e 2niv ' x 9. Then w 2 G V' A (M d+1 ) since for 
all k G {1, . . . ,d}, we have S Ak w 2 = e 2wiv < x+a ^®e = e 27riv a k e 27riv - x (g)e = l-e 27rivx (g>e = w 2 . 
Also, W2 G ^x>^(Rd+i)(p) because 

P (^? " ' ' cl) W2 = ( c «( 27riv ) + ci(2^v)6' + • • • + c fc (2^v)6( fc ) + o) e 2 ™ 

= 0-e 2mv - x = 0. 

Suppose it) patches up w\ := and w 2 , that is, w; G Q3yv(p) satisfies w|(_oo,o) = I (-00,0) = ® 
and w|( Ti +oo) = u>2|(r,+oo) 7^ 0- Proceeding as in the proof of the 'if part of Theorem 11.31 

w= <5v®rW, 

veA- 1 ^ 

for some T( v ) G Thus we have 

d f d \ ^ 

c (2mv)T + ci(2vriv) — T + • • • + c fe (27riv) - T = 0, 

and Tb^o) = 0, while Tb Tj+00 ) = 0|( T+oo ) 7^ 0, which is impossible, since Ck(2niv) ^= 0. 
This completes the proof of the fact that (1) =>■ (2). 

(2) =>■ (1): Let us suppose that Wi,w 2 G Q^^y (Rd+i) (p) • Then we have 



wi= Sv®t} v \ w 2 = ^v®T 2 (v) , 

for some T^\T 2 {v) gD'(K). Moreover, owing to the correspondence between D' A (M. d ) and the 
space of sequences s'(Z rf ) of at most polynomial growth, it follows that for each ip G T>(R), 
there exist M v > and a positive integer k<p, such that we have the estimates 

I ( T i V) , f) I < M^l + |n|)\ KT 2 (V) ,^)| < M„(l + |n|)\ 

for all n := Av G Z d . Let r > 0, and let G C°°(R) be such that 0(t) = 1 for all t < 0, 
6(t) = for all t > r/4 and < 6(t) < 1 for all i G R. Define G £>'(R) by 

T (v) ■- qt[ v) + 6{t - . 
Set u; G £>'(R d+1 ) to be w = ^ 5 V ® T (v) . Then for every ip G £>(K), we have 

vGA- 1 ^ 

^( v ) «M < l/flT_ (v) m\ \ -4- \ lft(r - -\ri v) 



|<TW^)| < 1(^,^)1 + |(0(r-.)rr^)l 

< M^(l + |n|) fe ^ + M, (r _. v (l + |n|) A 

< max{M $lp , M e{T _. )tp }{l + | n |) max ^' fc « 



>(t— )*> 



(t—)<p) 
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V 



and so w(<p) G V' A (R d ). Thus w G P^(M a!+1 ). Also, w G 55^ (M d+i)(p) because 

c (27Tiy)(5 v ® rW) + Cl (2™ y )^(<5 v ® T^) + • • • + c n (2vriy) (^)Vv ® T^) 
= co(2vriv)5 v ® T (v) + ci(27riv)5 v ® ^T (v) + • • • + c n (2mv)<5 v ® (^"t (v) 

= o, 

w = 0. 



and so 



c (2-Kiy)w + ci{2my)—w H h c n (27riy)(— ^ 

Finally, because T( v ) | (.^o) = zf v) | (-00,0) and ^ (v) l(r,+oo) = r 2 V) I(t,+oo)> ft follows that 
^l(-oo,o) = ^i|(-oo,o) and u>|( Tj+00 ) = io 2 |( Tj+00 ). Consequently, H(-oo,0) = ^i|(-t»,o) and 
w\(t,+oo) = ^2|(r,+oo) 5 showing that the behaviour is controllable. This completes the proof. 

□ 

Example 3.1 (Diffusion equation). Consider the diffusion equation 

d d 2 d 2 



dt dx\ ® x \! ' 

corresponding to to the polynomial p = r — £ 2 — ■ ■ ■ — £ j . Here the ^-content (p) is the ideal 
in C[£, . . . , £d] generated by the two polynomials 1 and — £f — — and so C^{p) is the full 
ring C[£, . . . , Consequently, its variety in C d is empty. Hence the behaviour %$T>' A {«. d + 1 )(p) 
is always autonomous, no matter what A is. Note that this is in striking contrast to what 
happens when we look at just distributional solutions: since 

deg(p(&, r)) = deg(r _(£ + ...+ £})) = 2 

Ik 

deg(p(0,...,0,r))=deg(r) = 1 

we have from Proposition 11.21 that i Sv'(R d + 1 )(p) * s n °t autonomous, and this outcome is 
physically unexpected. Indeed, if we imagine the case of diffusion of heat, in which case the 
w is the temperature, say along a metallic rod when d = 1, then zero temperature upto time 
t = should mean that the temperature stays zero in the future as well (since the above PDE 
describes the situation when no external heat is supplied). However, when one considers 
distributional solutions, one can have pathological solutions with a zero past that are nonzero 
in the future! But if we choose the physically "correct" solution space in this context, namely 
functions which at each time instant have a spatial profile belonging to L°°(M), then it can be 
shown that solutions that are zero in the past are also zero in the future, as expected. So the 
real reason for the nonautonomy when one considers solutions in V(M. d+1 ) is that there is no 
restriction on the spatial profiles of the solutions at each time instant, and wild growth (such 
as something which grows faster than e' x ' ) is allowed. However, with a periodic profile in 
the spatial direction, namely when the spatial profile is in D' A (U. d ), we know that the spatial 
profile is automatically tempered (see for example [5]), and as we have seen above, in this 
case the behaviour %$T>' A (R d + 1 )(p) ^ s autonomous, in conformity with our physical expectation. 

Similarly, since cq = 1 7^ 0, it follows from Theorem II .51 that (Rd+i)(p) is never control- 
lable. 
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Remark 3.2. We end with an open question: Is there an algebraic-geometric characteri- 
zation in terms of p of approximate controllability of Here, by approximate 
controllability, we mean the following. 

Let V' A (R d x (0,+oo)) := {w\ t>0 : w G V' A (R d+1 )}, endowed with the induced topol- 
ogy from P'(R d x (0, +00)). We call %$T>' A (R d + 1 ){p) approximately controllable if for every 
wi,W2 G and every neighbourhood N of in "D^(R d x (0, +00)), there is a 

w G 33-p/^ R d+i)(p) and a r > such that 

H(-oo,0) = w i|(-oo,o) and 
S (0,-t){w - W 2 )|( ,+oo) G N. 

Here, S( 0i - r ) denotes the translation operator in £>'(M d+1 ) corresponding to (0, — r) G 

When one considers just smooth solutions, that is, W = C°°(M. d+1 ), then one can define 
approximate controllability analogously to the above. (In particular, then the N in the 
definition is a neighbourhood in the appropriate topology of the Frechet space of C°° functions 
in the half-space t > 0.) Using the main result in [6], a characterization of approximate 
controllability in the smooth solution case was given in [9J Theorem 4.6]. 

It is not hard to show that for a nontrivial behaviour ©jy (r<j+i)(p), one has the following 
heirarchy of properties: 

controllability => approximate controllability -1 (autonomy). 

Thus in our search for the appropriate algebraic-geometric condition characterizing approx- 
imate controllability of 58^/ (R<*+i)(i?)> we expect an algebraic-geometric condition that lies 
between the two characterizations of controllability and approximate controllability given in 
this article in Theorems 11.31 and 11.51 
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